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ABSTRACT 

A significant fraction of novae exhibit a series of rebrightenings on the dechne branch of their hght 
curves. We use visual observations to study this phase in several well-observed novae. We find that 
these rebrightenings are isolated flare-like events on otherwise smooth light curves and we show that 
in most novae in our sample the time intervals between consecutive flares gradually increase as a 
geometric series; rebrightenings are equally spaced in logarithmic time. We also find a correlation 
between the rate of increase in the time between rebrightenings with the speed class of the nova in the 
sense that slower novae tend to increase the time intervals faster. We attribute these rebrightenings 
to instabilities in the envelope hydrogen burning and we mention other cases of such timing pattern 
in astronomy and natural sciences in general. 

Subject headings: instabilities - novae, cataclysmic variables - nuclear reactions, nucleosynthesis, abun- 
dances 



1. INTRODUCTION 



Decay light curves of explosive phenomena are often 
bumpy. Optical and X-ray light curves o f gamma-ray 
burst s are superimp osed with flare s (e.g. IStanek et al " 
2007tlDai et al.ll2007t l&jrrows et al.ir2005HO'Brien et al 



20061 ) suggesting late-time activity of the central engine. 
Similar features have been observed during outbursts 
of X -ray novae (e.g. ICallanan et"al] Il995l : IChen et all 
Il997f) and of an acc retion-driven millisecond pulsar 
(IWijnands et al.l l2001[ ). A fraction of classical novae 
also exhibit a series of rebrightenings that were observed 
for the flrst tim e during the 1901 outburst of GK Per 
(|Campbel]|ll903f l. These rebrightenings are the subject 
of this Letter. 

An outburst of a classical nova is a result of a ther- 
monuclear runaway in the matter accreted onto the sur- 
face of a white dwarf, whic h ejects an ex panding en- 
velope (for a review see e.g. IWarneijll995f ). Maximum 
visual brightness occurs when the pseudo-photosphere 
reaches its maximum radius. Then, the pseudophoto- 
sphere shrinks and heats while keeping the bolometric 
flux essentially constant. This results in a decrease in 
the visual flux as the bolometric luminosity is radiated 
at shorter wavelengths. After the ejection the hydrogen 
remaining on the white dwarf continues burning near the 
Eddington limit for several weeks to several years (e.g. 
Prialn ik 1986; Ogelman et "aIlll993HKrautter et al.lll996l : 
Ness et al.ll2n07ll2008[ ). 



Roughly a third of very fast or fast novae (|Warnejll995L 
p. 259) decline smoothly. The rest either show a series of 
rebrightenings or they create dust which brings the visual 
flux to a deep minimum. This time period, which usu- 
ally occurs 3-4 mag below visual maximum, is commonly 
called the "transition phase" because of spectroscopic 
changes toward a nebular spectrum. For simplicity, we 
will restrict the term "transition phase" throughout the 
Letter only to the rebrightenings. 
The underlying physical reason for the transition phase 
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is not yet clear, for a recent review of the proposed 
models, see ICsak et al.l ()2005() . They spectroscopically 
monitored the transition phase of V4745 Sgr and found 
that during the rebrightenings the continuum emission 
rose and that the spectral lines showed P-Cyg profiles, 
which suggests a temporary increase in mass-lo s s rate . 
On the basis of these observations, ICsak et all (|2005l ) 
suggested that the rebrightenings are envelope-hydrogen- 
burning instabilities. They also observed a gradual 
increase of the time intervals between rebrightenings, 
w hich was also not e d pre viously in GK Per and DK Lac 
bv lBianchini et all ()1992f ). 

In this Letter, we analyze visual light curves of no- 
vae that experience the transition phase. We show that 
these rebrightenings are indeed of eruptive nature and 
that they follow a specific pattern that relates to the 
overall rate of decline of the light curve of the nova. 
We conclude that the rebrightenings are likely caused 
by hydrogen-burning instabilities and mention other ex- 
amples of similar timing pattern in astronomy and in the 
theory of self-organizing critical systems. 

2. VISUAL LIGHT CURVES 

To study the short timescale variations during the tran- 
sition phase, we need well-sampled light curves obtained 
using a standard-enough detector so that we can directly 
compare different objects that erupted during the last 
~100 years. Due to the plethora of variable and often 
strong emission lines, a slight difference in the detec- 
tor response function could create substantial system- 
atic shift in the magnitude. We choose to use visual 
observations carried out by amateur astronomers. We 
obtained magnitude estimates from the archives of the 
AAVSO^ (Henden 200^ and AFOEV^ for V603 Aql, 
V1494 Aql, DK Lac, V2540 Oph, and V4745 Sgr. For 
V4745 Sgr, we supplement these data with visual obser- 
vations from the VSNET'^ archive and with the V band 
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Fig. 1. — Light curves of novae in our sample smoothed using moving averages with bin size of 1 day. The data were shifted so that the 
maximum occurred at tmax with brightness mag. V4745 S gr is plotted for refere nce in both panels. The right panel shows also X-ray 
light curve of accreting millisecond pulsar SAX J1808. 4-3658 llWiinands et al.ll200lh . 



measurements of the All Sky Automated Survey (AS AS; 
lPoimanski|[200l. For GK Per, we use the mean light 
curve of iCampbelll (fl90l constructed from visual ob- 
servations. We als o include V373 Set and its visual light 
curve published bv lRosinol (1X978') and supplement it with 
AAVSO and AFOEV observations. To our knowledge, 
the above mentioned objects represent the best-observed 
novae exhibiting sufficient number of rebrightenings to 
study their timing evolution. Note that the error in a 
single visual estimate i s usually assumed to be about 0.3 
mag (jKiss et al.l |1999() and we investigate events with 
amplitudes typically greater than 1 mag and defined by 
many observations of several observers. The rebrighten- 
ings we discuss here are robust. 

In Figure [U we show the visual magnitudes of novae in 
our sample as a function of time. To illustrate a nova 
light curve with a sm ooth decline, we include a very 
fast nova, V1500 Cyg (|Mattei|[l993[ ). We plot moving- 
averaged data with bin size of 1 day, and we align the 
light curves so that the visual magnitude at the time of 
maximum imax is mag. This procedure is potentially 
problematic in the case of V4745 Sgr and V2540 Oph, 
because the rebrightenings began shortly after discovery, 
unlike other novae. However, prediscovery ASAS obser- 
vations of V4745 Sgr rule out the possibility that the true 
visual maximum was missed. Given the overall similarity 
of light curves of both objects, we suggest that the same 
is true for V2540 Oph and we set tmax for both objects 
to be approximately the time of discovery. Importantly, 
an uncertainty of a couple of days in tmax does not affect 
any of our conclusions. 

Shortly after maximum, the visual fl ux of novae de- 
creas es as a power-law function of time (jHachisu fc Katol 
[2005h . which can be observed as a linear decline in our 
Figure [TJ Several hundred days after iniax7 the slope 



of the decline steepens, which is caused by a decrease 
in the wind mass-loss rate. The shell hydrogen burn- 
ing stops shortly afterward (|Hachisu fc Katoll2006l ). The 
rebrightenings start to occur approximately 20 days af- 
ter the maximum. At this time, fast novae like GK Per 
or V1494 Aql have faded by about 4 visual magnitudes 
while the slow novae V2540 Oph and V4745 Sgr are still 
close to their maximal visual brightness. 

To investigate the nature of rebrightenings more 
closely, we plot in Figure [5] a blow-up of well-covered 
rebrightenings for V4745 Sgr, DK Lac ,and GK Per. We 
see that the rise of brightness is much faster than the 
decline back to the base level. Further, for all novae the 
amplitude of the rebrightenings ranges from 1.0 to 2.0 
mag, without any obvious time dependence, although the 
duration of individual rebrightenings ranges from 2 days 
for GK Per to 10 days for V4745 Sgr. We note that even 
for GK Per, where the time interval between the max- 
ima is about 5 days, the individual rebrightenings are 
well separated by a period of essentially constant mag- 
nitude. It should be noted, however, that the data for 
V603 Aql and V1494 Aql do not show the rebrighten- 
ings so clearly separated. We attribute this to the lower 
quality of data available for these stars and to the aver- 
aging that smears out details in the light curves. Indeed, 
the phase of constant brightness stands out more clearly 
in unbinned data. Given this evidence, we consider the 
transition phase rebrightenings to be appropriately de- 
scribed as flares on top of a continuous underlying light 
curve. 

3. LONG-TERM BEHAVIOR OF THE FLARES 

Looking again at Figure [T] and particularly at the light 
curve of V4745 Sgr (purple line) we see that the flares 
have approximately constant separation in logarithmic 
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Fig. 2. — Details of individual rebrightenings. Data for GK Per 
are plotted with green triangles, for DK Lac with black stars, 
and visual data (purple open circles) for V4745 Sgr are plotted 
along with ASAS V-band measurements (purple squares). The 
light curves were aligned so that the rebrightenings start at tstart 
and arbitrary shifts in magnitude were applied for illustration pur- 
poses. Purpose of the solid lines is to guide the eye when reading 
the plot. 

time. This means that the intervals between consecutive 
m axima get graduall y longer, as was noted for this nova 
bv ICsak et atl (|2005l) . V2540 Oph and DK Lac exhibit 
apparently the same behavior. To investigate further this 
gradual increase of separation between flares, we have de- 
termined times of maxima U, i = 1, 2, 3, . . ., by fitting a 
low-order polynomial to the data with the flare and read- 
ing its maximum position. We estimate the error of ti to 
be 10% of the time interval used for polynomial fitting. 
This translates to less than 1 day for most of the flares. 
We experimented with fitting data averaged with differ- 
ent bin sizes and with unbinned data and conclude that 
the differences in maxima timings U are smaller than 
the assumed error. We c ompared our maxima timings 
of DK Lac with those o f iBochn fcck (1951) and found 
good agreement. Data of iBochiifcek (,195L) contain tim- 
ings of several maxima that occurred prior to those that 
we were able to derive from our data. We add these 
timings to ours but we note that our conclusions below 
are not affected in any way by their inclusion. We also 
point out that in addition to the seven flares observed by 
ICsak et al.l (|2005[) in V4745 Sgr the ASAS data contain 
two additional flares that occurred at JD 2453165 and 
JD 2453185 {t - t^ax = 423 and 443 days) which have 
lower amplitude and shorter duration than the preced- 
ing flares. Two single brighter data points indicate that 
there might have also been other flares at JD 2453085 
and JD 2453420 (t - t^ax = 343 and 678 days) but we 
cannot show that they are maxima. 

In the left panel of Figure [H we show the time interval 
between two consecutive flares ti — ti-i as a function 
of time elapsed since i,nax for V4745 Sgr, DK Lac, and 



GK Per. We see that the timings for V4745 Sgr and 
DK Lac follow a power law 

\og{ti - U-i) = a + b\og{ti - imax) (1) 

with similar values of the power-law index b but with 
different normalization a. A least-squares fit to the data 
yields b = 0.88 ± 0.04 for DK Lac and b = 0.79 ± 0.04 
for V4745 Sgr, assuming a 2-day error of maxima tim- 
ings. However, maxima timings for these stars can be 
well fitted with b being held fixed at 6 = 1 as is shown in 
Figure O This value of b has a convenient property that 
the ratio of the times of two consecutive flare maxima 
is constant: {U - tmax)/(ij-i - ^max) = 1/(1 - 10°). In 
other words, the times of the flares ti — tmax increase as 
a geometrical series. 

The situation is more complicated for GK Per. Ini- 
tially, th e time interval betwee n the flares increases, as 
noted bv lBianchini et all ()1992D . and the increase agrees 
with Equation ([1]) with b = 1 (as is shown by a fit to the 
first eight points in the left panel of Figure [3|) . After- 
ward, the flare separation remains steady with a possible 
decrease at ti — tmax ~ 100 days. While this behavior 
is not seen in DK Lac, we point out that the two latest 
flares of V4745 Sgr are separated by only about 20 days 
and suggest that a late decrease of flare separation can 
occur also in novae that otherwise behave according to 
Equation ([T]). 

Time intervals between consecutive flares in V4745 Sgr 
grow faster than in DK Lac, the normalization a of the 
former is higher than that of the latter. Flare intervals for 
GK Per are consistently lower than for both DK Lac and 
V4745 Sgr at the same t — tmax- Is there some correlation 
between this behavior and other parameters of the novae? 
Looking at Figure [U we find that the visual magnitude 
of GK Per declined faster than that of DK Lac which is 
in turn slower than V4745 Sgr. We propose there is a 
correlation between the normalization a and the rate of 
decline of the nova. 

To check this assumption further, we plot in the right 
panel of Figure [3] time intervals between consecutive 
flares for the remaining novae in our sample plus those 
flare timi ngs of V 3 73 Set determined from the light curve 
given by iRosinol (|1978D that are supported by obser- 
vations. We see that V2540 Oph with basically the 
same light curve follows almost the same power law as 
V4745 Sgr, albeit with a lower total count of flares. Sim- 
ilarly, flare times of V373 Set fall between DK Lac and 
V4745 Sgr, which is true also for the light curve. This 
correlation does not entirely hold for V1494 Aql and 
V603 Aql. V603 Aql is usually regarded as a very fast 
nova similar to GK Per and V1494 Aql. However, while 
its visual light curve declines initially in the same man- 
ner as GK Per, starting from t — tmax ~ 50 days it de- 
clines less steeply and in the final stages of the light curve 
V603 Aql is slower than GK Per by a factor of about 2. 
On the other hand, the light curve of V1494 Aql is vir- 
tually the same as of GK Per, but the flare separations 
are significantly larger and until t — i„iax ~ 90 days they 
increase according to Equation ([T]). Interestingly, the to- 
tal duration of the transition phase is the same in both 
objects. 

To put our correlation analysis on more quantitative 
grounds, we have obtained the rate of the magnitude 
decline by determining the time-stretching factor p of 
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Fig. 3. — Time intervals between two consecutive flares ti — ti—\ as a function of time elapsed since the optical maximum t j — imax • Left 
panel: data for V4745 Sgr (purple circles), DK Lac (black stars), and GK Per (green triangles pointing up). Lines mark fit of Equation |(T]| 
with 6 = 1 to the all points of DK Lac, first six points of V4745 Sgr and first eight points of GK Per. Right panel: the same as in the 
left panel but for V2540 Oph (orange squares), V1494 Aql (blue triangles), V603 Aql (open red diamonds) and V373 Set (black half-filled 
squares). Overplottcd are fits of V4745 Sgr, DK Lac and GK Per from the left panel. 

the light curves of our novae with respect to the hght 
curve of V1500 Cyg, assuming that hgh t curves of novae 
are homologous (jHachisu fc Katoll2006f ) . We also derived 
normalizations a while holding 6=1 fixed for all our no- 
vae. Spearman's rank correlation coefficient of quantities 
a and p is 0.75 with a 6% chance of random correlation. 
We conclude that the time intervals between the flares 
generally increase faster for slower novae, but there ex- 
ists other parameter or parameters that have also effect 
as is shown by differences in V1494 Aql and GK Per. 

Finally, we note that while the time interval between 
the flares increases, the duration of the individual flares 
does not increase in similar manner. For example, in 
GK Per the duration varies randomly by about a factor of 
2 with a mean of ~3 days and in DK Lac the duration of 
flares ranges from 4 to 9 days. This can be observed also 
in Figure [T] where the flares are becoming increasingly 
narrower due to the nonlinear abscissa. 

4. DISCUSSION 

During the shell hydrogen burning, when the transi- 
tion phase occurs, the white dwarf is radiating at its 
Eddington luminosity (Prialnik 1986). An isotropic en- 
ergy release of 10'''^-10 ergs is required to brighten the 
nova by an additional factor of ~5 for about 5 days, as- 
suming the rebrightening affects between 3% and 100% 
of the bolometric luminosity. This can be accomplished 
either by burning of 10~^-10~* of hydrogen, by ac- 
creting from 10~*-10~^ M0 of matter onto the white 
dwarf, or by tapping onto the rotat ional energy of the 
white dwarf (jHachisu fc Katol [2009D . We do not con- 
sider accretion likely because the white dwarf powers a 
conti nuum-driven wind during the shell hydrogen burn- 
ing (Prialniklll986D . Rotational energy extraction would 
happen most likely by employing strong magnetic fields 



embedded in the expanding envelope and rotating differ- 
entially with respect to the white dwarf. When the enve- 
lope density decreases enough, the magnetic field restores 
synchronous rotation driving strong ma gnetic activities. 
This process has been suggested by iHachisu fc Katol 
(|2009l) as an explanation for the secondary maxima of 
several recent novae. However, a strong magnetic field 
is not a sufficient condition for the rebrightenings to oc- 
cur: [Schmidt & Stockman (1991) found polar strength 
magnetic fields in a smoothly decaying nova V1500 Cyg. 

IShavivl (|2001h suggested that the transition phase in 
novae is caused by sta gnating super-Eddington winds. 
IVan Marie et all (|2009D performed numerical simulations 
of porosity-moderated photon-tiring winds and found 
strong quasi-periodic luminosity variations caused by the 
fact that a changing fraction of the photon energy is 
being used to drive the wind. They also found com- 
plex patterns of outflowing and infalling matter. We 
found that the transition phase flares are generally iso- 
lated events well separated by intervals of relatively con- 
stant or slowly de clining visual brightness. Furthermore, 
ICsak et al] ()2005() analyzed optical spectra covering the 
transition phase of V4745 Sgr and did not find any spec- 
troscopic evidence of infalling matter. They concluded 
that the flares are caused by repetitive instabilities in 
shell hydrogen burning on the white dwarf. Our esti- 
mate of 10~^-10~® M0 of hydrogen burned is 1-2 or- 
ders of magnitude lower than the minimum envelope 
mas s necessary to sustain equilibrium hydrogen burn- 
ing (|Sala fc Hernanj|2005[ ). A late-time increase in the 
envelope hydrogen burning and the associated mass ejec- 
tion has also been invoked to ex plain a secondary peak in 
V2362 Cyg ()Lvnch et al.l 120081) . Therefore, we conclude 
that the transition phase flares are likely caused by sud- 
den increases in otherwise steady hydrogen burning in 
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the white dwarf envelope. 

By analysis of flare timings, we have shown that the 
intervals between consecutive flares increase as a power- 
law function of time. However, the length of the flares 
remains essentially constant during the transition phase 
suggesting that only the trigger of the flare is affected. 
The power-law exponent is generally the same for all no- 
vae, but the normalization is larger for novae with smaller 
rates of decline. The primary parameter that controls the 
rate of decline is the mass of the white dwarf: faster novae 
occur on more massive white dwarfs (,Prialnik fc Kovet j 
II995I ) . It is then natural to expect that with a more mas- 
sive white dwarf the intervals between hydrogen-burning 
flashes will be shorter and also that the flare will last 
a shorter amount of the time. However, as we have 
demonstrated on V1494 Aql and GK Per, which have es- 
sentially the same light curve but substantially different 
intervals between flares, there is one or more parameters 
that also influence the properties of the flares. We sus- 
pect that chemical composition of the white dwarf and 
its hydrogen-burning envelope, as well as its mass or the 
mass accretion rate might play a role. 

The gradual increase of intervals between individual 
rebrightenings is an unusual phenomenon in astronomy. 
We have been able to find only one other example of 
such behavior: the X-ray and optical light curve of the 
accretion-driven millisecon d pulsar SAX J1808. 4-3658 
during its 2000 outburst (jWiinands et al.ll2001l ). We 
plot the X-ray light curve in the right panel of Figure [H 
However, during the 2002 outburst the flares were more 
regular, with no apparent increase in time separation 
(|Wiinandsll200l . 

In other areas of physics the phenomenon of logarith- 
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mic periodicity occurs as a correction to a power-law 
behavior a nd suggests pre sence of discrete scale invari- 
ance (DSI, ISornettd [200I p. 276). A well-known ex- 
ample of DSI is the triadic Cantor set which is geo- 
metrically identical to itsel f only und er magnifications 
equal to powers of 3 (.Sornettd Il998f l. Among exam- 
ples of DSI and logarithmic periodicity observed in na- 
ture we specifically mention the rate of acoustic emis- 
sions as a f unction of time before t he rupture of pres- 
sure vessels (jZhou &: Sornettdl2Q02 a^ , and ground- water 
chemical anomalies pr e ceding the 1995 Kobe ear thquake 
(jJohansen et all 120001 : IZhou &: Sornettd l2002H l. The- 
oretical models, which exhibit spontaneous emergence 
of DSI, include earthq uake aftershock simu lations based 
on cellular automata ( Lee fc Sorn ette 200 0|) and rupture 
studies with a thermal fuse model (Joha nsen fc Sornettd 
|1998[ ). Log-periodic echoing was also found during a crit- 
ical r elativistic gravitational collapse of a rn assless scalar 
field (|ChoDtuiklll99l IChoptuik et al.l[200l . Finally, we 
point out that there is no si ngle underlying mechanism 
for the DSI in general (|Sornette..2004 . p. 160). 
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